Next generation extended Lagrangian first principles molecular dynamics.
نویسنده
چکیده
Extended Lagrangian Born-Oppenheimer molecular dynamics [A. M. N. Niklasson, Phys. Rev. Lett. 100, 123004 (2008)] is formulated for general Hohenberg-Kohn density-functional theory and compared with the extended Lagrangian framework of first principles molecular dynamics by Car and Parrinello [Phys. Rev. Lett. 55, 2471 (1985)]. It is shown how extended Lagrangian Born-Oppenheimer molecular dynamics overcomes several shortcomings of regular, direct Born-Oppenheimer molecular dynamics, while improving or maintaining important features of Car-Parrinello simulations. The accuracy of the electronic degrees of freedom in extended Lagrangian Born-Oppenheimer molecular dynamics, with respect to the exact Born-Oppenheimer solution, is of second-order in the size of the integration time step and of fourth order in the potential energy surface. Improved stability over recent formulations of extended Lagrangian Born-Oppenheimer molecular dynamics is achieved by generalizing the theory to finite temperature ensembles, using fractional occupation numbers in the calculation of the inner-product kernel of the extended harmonic oscillator that appears as a preconditioner in the electronic equations of motion. Material systems that normally exhibit slow self-consistent field convergence can be simulated using integration time steps of the same order as in direct Born-Oppenheimer molecular dynamics, but without the requirement of an iterative, non-linear electronic ground-state optimization prior to the force evaluations and without a systematic drift in the total energy. In combination with proposed low-rank and on the fly updates of the kernel, this formulation provides an efficient and general framework for quantum-based Born-Oppenheimer molecular dynamics simulations.
منابع مشابه
Extended Lagrangian free energy molecular dynamics.
Extended free energy Lagrangians are proposed for first principles molecular dynamics simulations at finite electronic temperatures for plane-wave pseudopotential and local orbital density matrix-based calculations. Thanks to the extended Lagrangian description, the electronic degrees of freedom can be integrated by stable geometric schemes that conserve the free energy. For the local orbital r...
متن کاملEffects of asymmetric stiffness on parametric instabilities of rotor
This work deals with effects of asymmetric stiffness on the dynamic behaviour of the rotor system. The analysis is presented through an extended Lagrangian Hamiltonian mechanics on the asymmetric rotor system, where symmetries are broken in terms of the rotor stiffness. The complete dynamics of asymmetries of rotor system is investigated with a case study. In this work, a mathematical model is ...
متن کاملStable and Efficient Linear Scaling First-Principles Molecular Dynamics for 10000+ Atoms.
The recent progress of linear-scaling or O(N) methods in density functional theory (DFT) is remarkable. Given this, we might expect that first-principles molecular dynamics (FPMD) simulations based on DFT could treat more realistic and complex systems using the O(N) technique. However, very few examples of O(N) FPMD simulations exist to date, and information on the accuracy and reliability of t...
متن کاملPerformance of extended Lagrangian schemes for molecular dynamics simulations with classical polarizable force fields and density functional theory.
Iterative energy minimization with the aim of achieving self-consistency is a common feature of Born-Oppenheimer molecular dynamics (BOMD) and classical molecular dynamics with polarizable force fields. In the former, the electronic degrees of freedom are optimized, while the latter often involves an iterative determination of induced point dipoles. The computational effort of the self-consiste...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- The Journal of chemical physics
دوره 147 5 شماره
صفحات -
تاریخ انتشار 2017